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In recent years, there has been an intense interest in understanding the microscopic mechanism of 
thermally induced magnetization switching driven by a femtosecond laser pulse. Most of the effort 
has been dedicated to periodic crystalline structures while the amorphous counterparts have been 
less studied. By using a multiscale approach, i.e. first-principles density functional theory combined 
with atomistic spin dynamics, we report here on the very intricate structural and magnetic nature 
of amorphous Gd-Fe alloys for a wide range of Gd and Fe atomic concentrations at the nanoscale 
level. Both structural and dynamical properties of Gd-Fe alloys reported in this work are in good 
agreement with previous experiments. We calculated the dynamic behavior of homogeneous and 
inhomogeneous amorphous Gd-Fe alloys and their response under the influence of a femtosecond 
laser pulse. In the homogeneous sample, the Fe sublattice switches its magnetization before the 
Gd one. However, the temporal sequence of the switching of the two sublattices is reversed in 
the inhomogeneous sample. We propose a possible explanation based on a mechanism driven by 
a combination of the Dzyaloshiskii-Moriya interaction and exchange frustration, modeled by an 
antiferromagnetic second-neighbour exchange interaction between Gd atoms in the Gd-rich region. 

We also report on the influence of laser fluence and damping effects in the all-thermal switching. 


I. INTRODUCTION 

Switching the sublattice magnetization directions of 
amorphous Gd-Fe alloys [1] (doped with small amounts 
of Co) by intense femtosecond laser pulses has gener¬ 
ated significant interest both experimentally and the¬ 
oretically. Amorphous Gd-Fe alloys are ferrimagnetic, 
with a strong antiferromagnetic (AFM) coupling between 
the rare-earth and transition metal moments, a coupling 
which has its explanation in the hybridization of the 
5d and 3d orbitals of the constituting elements [2]. In 
Ref. [1], it was found that an optical excitation caused 
the net magnetization of both the Gd and Fe sublat¬ 
tices to rapidly collapse. However, the time scales of the 
dynamics of the two sublattices were found to be quite 
different: the net magnetic moment of the Fe sublattice 
was found to vanish after 0.4 ps and then for a short pe¬ 
riod of time, up to 2 ps, be parallel to the Gd moment. 
The Gd sublattice, which initially is antiferromagnetic to 
Fe, vanished after 2 ps, after which it reversed to have its 
magnetization opposite to that of Fe, hence completing 
the reversal process (see Fig. 3 of Ref. [1]). The interest 
of these results obviously have great potential for techno¬ 
logical applications, since they open up for possibilities 
to store information in a magnetic medium without ap¬ 
plying an external magnetic or electric field. In fact, the 
experimental results reported in Ref. [1] follow intense 
investigations of magnetization dynamics, which started 
in the mid ’90s [3-8] . 

Different theoretical models have attempted to explain 
these results. For instance, in Ref. [9] it was argued that 


two time and temperature domains were relevant, where 
the spin-relaxation was driven first by a relativistic con¬ 
tribution whereas subsequently relaxation was argued to 
be governed by an exchange origin. A different explana¬ 
tion was provided in Ref. [10] where the coupling between 
Gd and Fe dominated magnon modes were identified as 
the most important aspect of the complex switching be¬ 
havior of the Gd-Fe system. It should also be mentioned 
that in the experimental investigation of Ref. [11], it was 
speculated that angular momentum was transferred be¬ 
tween the different sublattices via spin-currents, and this 
was identified as the most important aspect of the mag¬ 
netization dynamics of amorphous Gd-Fe alloys. 

Although the main experimental findings of Ref. [1] 
have been repeated in subsequent experiments, there are 
details in a more recent work that have so far not been ad¬ 
dressed satisfactorily by theory. For instance, in Ref. [11] 
several hitherto unexplained experimental facts were re¬ 
ported. Moreover, in the samples measured by Graves 
et al. [11] concentration profiles were detected, with Gd 
rich/Fe poor regions and Gd poor/Fe rich regions, in the 
same sample. Surprisingly, it was found that for the Gd 
rich regions, the Gd moment has a different dynamical 
response compared to that of the Gd poor regions. This 
amounted to situations where in the Gd rich regions, the 
Gd moment reversed before the Fe moment, in contrast 
to the result reported for the average magnetization of 
Gd or Fe sublattices, reported in Ref. [1]. Hence, in 
amorphous Gd-Fe alloys, it seems that sometimes the Fe 
moment reverses before the Gd moment, and sometimes 
Gd switches before Fe, depending on the concentration 
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of Gd and Fe in local regions of the sample. 

None of the theories presented so far has addressed the 
role of the amorphous structure and the chemical inho¬ 
mogeneity of the Gd-Fe system and how this influences 
the ultrafast switching behavior. In this work, we present 
a multi-scale approach to address this problem, where 
we coupled first principles electronics structure theory to 
atomistic spin dynamics simulations [12]. After intro¬ 
ducing our methodology, we substantiate our approach 
by comparing equilibrium magnetization curves with ex¬ 
periment for a wide range of concentrations, elucidating 
the crucial role of the amorphous atomic arrangement on 
the magnetism. Subsequently, we demonstrate that the 
results for homogenous samples are in agreement with 
previous theoretical analysis reported in Ref. [9]. Finally 
we turn to inhomogenous samples and demonstrate that 
the switching is crucially affected by the chemical inho¬ 
mogeneity and the non-collinearity of the spins in the 
rare-earth sublattice. In the appendix, the methods are 
explained in more detail and an analysis of the role of 
the damping is also provided. 


II. METHOD 

A. Details of the simulation of structural 
properties 

First-principles spin polarized calculations were per¬ 
formed by means of the density functional theory [13, 14] 
and projector augmented wave [15, 16] method as im¬ 
plemented in the Vienna ab initio simulation package 
(VASP) [17-19]. The exchange-correlation potential 
was treated using the generalized gradient approxima¬ 
tion with the Perdew, Burke, and Ernzerhof functional 
[20], including the valence states for Gd 

and 3(f4:S^ for Fe. The LDA+U method [21] was applied 
to Gd with Ueff=7 eV and J=1 eV. 

The amorphous structures were generated by means 
of the stochastic quenching (SQ) method [22, 23], as de¬ 
scribed in Ref. [24] . This method is based on the single¬ 
random-valley approximation in vibration-transit (VT) 
theory [25, 26]. The SQ method was demonstrated to 
provide reliable atomic coordinates of amorphous mate¬ 
rials [24]. In the initial structures, 200 atoms were both 
spatially and chemically randomly distributed in a cubic 
unit cell with periodic boundary conditions and a den¬ 
sity of g=7.S7 g/cm^ for Gdo. 24 Feo. 76 , g=7.SS g/cm^ for 
Gdo. 5 oFeo. 50 , and p=7.89 g/cm^ for Gdo. 76 Feo. 24 - The 
atomic positions were then relaxed until the force on ev¬ 
ery atom was negligible. The calculations were performed 
using the F point. 

The kinetic energy cutoff of 550 eV together with 
Methfessel-Paxton band smearing [27] of cr = 0.2 eV were 
used for electronic structure calculations. The atomic 
charges were determined from Bader analysis [28-30]. 


B. Details of the atomistic spin-dynamics 
simulations 

In our simulations, we combined the two-temperature 
(2T) model [31] with the atomistic spin dynamics (ASD) 
in the UppASD code [12] using the Landau-Lifshitz- 
Gilbert (LLG) equation. Model exchange parameters 
were used for all simulations. At a finite temperature, 
the temporal evolution of individual atomic moments in 
an effective field is governed by Langevin dynamics, 

= -7mjX[Bj + bj]-7—miX(miX[Bj + bj]), (1) 

at m 

where 7 is the gyromagnetic ratio, a represents the di¬ 
mensionless phenomenological Gilbert damping constant 
and mi stands for an individual atomic moment on site i. 
The “effective” magnetic field is represented by Bi while 
hi is a time evolved stochastic magnetic field, which de¬ 
pends on the electron temperature from the 2T model. 
After applying a femtosecond laser pulse on the sam¬ 
ples, the electron temperature increase from the initial 
temperature Tq to a peak temperature in less than 50 
fs. Then, the electron temperature slowly cools down in 
about 5 • 10^ fs since the heat of the electron system is 
transferred to the phonon bath via electron-phonon in¬ 
teractions [31]. With this method, details of all thermal 
switching are investigated in detail, and reported upon 
below. 


III. RESULTS AND DISCUSSION 
A. Ab-initio theory and structural properties 

In this section we provide structural information of 
Gda,Fei_a, {x = 0.24, 0.50, 0.76) magnetic alloys based on 
ab initio theory. We note here that the first-principles 
calculations result in a metallic character of these amor¬ 
phous materials, in agreement with experimental obser¬ 
vations. The electronic properties of Gd-Fe systems are 
described in more details in the Appendix A. 

The local atomic environment in amorphous Gda,Fei_a, 
can be analyzed by using radial distribution functions 
(RDF) calculated for different atomic pairs (see Fig. 1). 
From RDFs of the SQ-generated structures, we find 
short-range order up to 8 A for Gd-Gd and Gd-Fe, and 
up to 6 A for Fe-Fe atomic pairs. Gd-Gd, Gd-Fe and 
Fe-Fe bond lengths, extracted from RDFs, are shown in 
Table I along with the bond lengths in selected reference 
systems. The theoretical Gd-Gd bond length in amor¬ 
phous Gda,Fei_a, is found to be shorter, and therefore the 
bonds are stronger than in hexagonal close-packed Gd. 
At the same time, the Gd-Gd bond length is longer than 
in crystalline compounds consisting of Fe and Gd, such 
as cubic GdFe 2 , and trigonal GdFes, suggesting a weaker 
bonding in the amorphous matrix. We find a favorable 
agreement between theoretical Gd-Gd bond distance in 
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the Fe-rich {x = 0.24) and equiatomic {x = 0.50) amor¬ 
phous systems compared to the experimental values for 
melt-quenched amorphous Gdo. 22 Feo .78 and Gdo. 56 Feo. 44 , 
respectively (see Table I). The theoretical Gd-Fe bond 
length in amorphous Gda,Fei_a, is very close to that in 
crystalline GdFe 2 and GdFcs. We find the bond distance 
between Fe atoms in Gda,Fei_a, to be shorter than in bee 
Fe and cubic GdFe 2 , but at the same time larger than 
in trigonal GdFes. We also find a remarkable agree¬ 
ment between theoretical bond lengths between differ¬ 
ent pairs of atoms, such as Gd-Gd, Gd-Fe and Fe-Fe, in 
Gda,Fei_a, {x = 0.24,0.50) obtained by SQ simulations 
and the experimental ones for quench-melted Gda,Fei_a, 
{x = 0.22,0.56) with similar stoichiometry. This illus¬ 
trates the efficiency and accuracy of the SQ method to 
describe the structural properties of amorphous materi¬ 
als. 
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FIG. 1. (Color online) Radial distribution function in 
GdxFei-cc with three different stoichiometries. The dashed 
lines show data calculated with DFT, while the black solid 
line represents data provided by molecular dynamics calcula¬ 
tions. 


Next, we analyze the local environment in amorphous 
Gda,Fei_a,, as it is represented by the average coordi¬ 
nation numbers (Table II). With increasing Gd concen¬ 
tration, the Gd-Gd average coordination number, as ex- 


TABLE 1. Bond lengths (A) in amorphous Gd^Fei-a, system. 
The bond distances in selected crystalline and amorphous sys¬ 
tems are listed for comparison. 


System 


Gd-Gd Gd-Ee Ee-Ee 

Gdo. 24 Eeo .76 


3.47 

3.02 

2.50 

Gdo.5oEeo.5o 


3.55 

2.98 

2.43 

Gdo.76Eeo.24 


3.52 

2.97 

2.44 

hep Gd [34] 


3.57 

- 

- 

bee Ee [37] 


- 

- 

2.54 

GdEes [35] 


3.21 

2.97 

2.37 

GdEe2 [36] 


3.22 

3.08 

2.63 

am-Gdo. 22 Eeo .78 

[32] 

3.47 

3.11 

2.57 

am-Gdo.56Eeo.44 

[32] 

3.54 

2.95 

2.51 


pected, increases from 4.6 for x = 0.24 to 10.9 for 
X = 0.76. Similarly, the Gd-Fe coordination number in¬ 
creases with the number of Gd atoms. In case of the Gd- 
Fe atomic pair, the coordination number is almost five 
times smaller in the Fe-rich amorphous matrix compared 
to the Gd-rich one. Also, while the Fe concentration de¬ 
creases within the series (from x = 0.24 to x = 0.76), so 
does the Fe-Fe coordination number. However, the total 
coordination numbers on both Gd and Fe coordination 
shells overall change similarly, i.e. decrease through the 
series. The reduction in the total coordination number 
both for Gd and Fe atoms can be referred to the change 
from a more close packed structure (Fe-rich system) to a 
more open one (Gd-rich system). 


TABLE II. Average coordination numbers for amorphous 
Gda,Eei-cc system. Goordination numbers in reference sys¬ 
tems are listed for comparison in the lower part of the table. 
In GdEes the coordination for Gd and Ee atoms with different 
site symmetries is different. Therefore, we show coordination 
for all inequivalent positions (specified in parentheses). 


System 


Gd-Gd 

Gd-Ee 

Ee-Gd 

Ee-Ee 

Gdo. 24 Eeo .76 


4.6 

11.0 

3.5 

7.6 

Gdo.5oEeo.5o 


9.1 

5.7 

5.7 

4.5 

Gdo.76Eeo.24 


10.9 

2.4 

7.5 

2.0 

GdFea [38] 


4 

12 

6 

6 

GdFes [35] 


(3 a) 2 
(6c) 1-3 

(3a) 6-12 
(6c) 3-6 

(36) 6 (36) 6 

(6c) 3 (6c) 3 

(186) 1-2 (186) 1-2 

am-Gdo. 22 Feo .78 

[32] 

3.0 

8.8 

2.5 

7.9 

am-Gdo.56Feo.44 

[32] 

7.5 

3.2 

4.2 

3.0 


B. Generation of amorphous samples using 
molecular dynamics 

As commented in Sec. Ill A, we optimised the struc¬ 
tures of Gda,Fei_a, (x=0.24, 0.50, 0.75) magnetic alloys 
by means of ab initio methods considering a supercell 
of 200 atoms, but the lack of crystal periodicity in amor¬ 
phous structures led us to consider even bigger supercells 
just to be sure that the results are reliable and the physics 
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of the amorphous structure was fully captured. The size 
of the uew supercells is beyoud the limits of the preseut 
state of the art of the computatioual resources usiug a 
DFT methodology. lu order to deal with bigger super¬ 
cells, we employed a molecular dyuamics approach, so 
that the dyuamics of atomistic Fe aud Gd spius showu iu 
upcomiug sectious have beeu performed usiug as iuput 
parameters the structural data provided by the molec¬ 
ular dyuamics method. Cousequeutly, we adapted a 
two step procedure. luitially, we coustructed a cubic 
uuit cell (1600 atoms) with Gd aud Fe atoms usiug a 
deuse-raudom-packiug-of-hard-spheres (DRPHS) model 
aud usiug as iuput the lattice parameter ( 29.5 A) pro¬ 

vided by geometry optimizatiou iu VASP calculatious 
for a cell of 200 atoms [39]. Iu the secoud step, the 
forces ou atoms were miuimized by usiug au adaptiug 
Morse pair poteutial iu Large-scale Atomic/Molecular 
Massively Parallel Simulator (LAMMPS) code [40]. Dur- 
iug the optimizatiou process, experimeutal boud leugths 
betweeu the species aud poteutial depths were takeu from 
Ref. [41]. The relaxed molecular dyuamics (MD) amor¬ 
phous samples were aualyzed via their RDF maiu peak 
positious aud uearest-ueighbour distributious. We illus¬ 
trate iu Fig. 1 the radial distributiou calculated by DFT 
aud MD simulatious for Gdo. 24 Feo. 76 - It may be observed 
that the agreemeut is rather good betweeu the two sets 
of theoretical values. 


C. Magnetization-dynamics and all-thermal 
switching 

1. Curie and compensation temperatures 

To further exteud the applicability of our method¬ 
ology, we performed ASD simulatious usiug the Up- 
pASD method ou a cell coutaiuiug 1600 atomic spius 
with periodic bouudary couditious. We take a ferri- 
maguetic amorphous GdFe model system based ou the 
structural parameters provided by both molecular dy¬ 
uamics aud DFT calculatious. The microscopic model 
exchauge parameters were takeu from Ref. [42]. We 
used the bulk exchauge values for ueighbouriug TM aud 
RE ious (Ji7e-Fe=0.8 mRy, JGd-Gd=^‘^^ mRy) because 
they provide the correct Gurie temperature for the re¬ 
spective sublattices. The value of the iutersublattice ex¬ 
chauge coupliug {Jcd-Fe ="0.25 mRy) was choseu to fit 
the temperature depeudeuce of the saturatiou magueti- 
satiou of both Fe aud Gd sublattices with results of x-ray 
maguetic circular dichroism measuremeuts of static mag- 
uetisatiou reported iu Ref. [42]. For the maguetic mo- 
meuts we take the bulk values, i. e. 7.6 fiB aud 2.1 {jlb 
for Gd aud Fe, respectively. No exterual field was applied 
iu the simulatious aud the auisotropy of the Gd aud Fe 
sublattices was ueglected iu the hamiltouiau. The model 
exchauge parameters are capable of reproduciug the key 
static maguetic properties, especially the Curie temper¬ 
atures aud maguetic compeusatiou poiuts. To illustrate 


this fact we show iu Fig. 2 both the calculated compeusa¬ 
tiou aud Curie temperatures for differeut Gd couceutra- 
tious, iu the rauge 20 at.% to 30 at.%. The Curie tem¬ 
peratures have beeu calculated usiug a fiuite size scaliug 
aualysis as described iu Ref. [43] (see Appeudix B). We 
fouud a very good agreemeut with the reported experi¬ 
meutal results. Moreover, we observe a geueral treud for 
the Tc to decrease as the Gd couceutratiou iucreases. 
We attribute this maguetic softeuiug to the additiou of 
more Gd-Gd uearest ueighbours, which have a smaller 
exchauge coupliug as compared to the Fe-Fe iuteractiou. 
But still, siuce we have au amorphous structure aud a 
supercell with limited uumber of atoms, this couditiou is 
uot fulfilled for every Gd couceutratiou aud this is why 
at some couceutratious, the Tc cau still iucrease slightly, 
as is the case for a couceutratiou of Gd of about 25 at.%. 
Iu the thermodyuamic limit, these smaller fluctuatious 
of the Tc are expected to vauish. It may be seeu from 
Fig. 2 that the simulatious reproduce with good accuracy 
the measured compeusatiou temperatures as well. Both 
the measured treud aud the absolute values of the com¬ 
peusatiou temperature of these alloys are reproduced by 
theory, where the most uoticeable feature is the iucrease 
of the compeusatiou temperature with iucreasiug Gd cou¬ 
ceutratiou. The reasou for this treud is a competitiou of 
maguetic sublattices with autiparallel coupliug. Too few 
Gd atoms result iu a Gd sublattice with a uet mague- 
tizatiou that is smaller thau that of the Fe sublattice, 
already at low temperatures, aud there is uo compeusa¬ 
tiou poiut. However, with iucreasiug Gd couceutratiou, 
the uet maguetic momeut of this sublattice is larger thau 
that of the Fe sublattice, at low temperatures. Siuce 
the Fe exchauge is strouger thau the Gd exchauge, the 
Gd sublattice maguetizatiou decays faster with temper¬ 
ature compared to the Fe sublattice maguetizatiou, aud 
at the compeusatiou temperature, they have equal size 
aud opposite directiou. Iucreasiug the Gd couceutratiou 
makes the maguetizatiou of this sublattice strouger rel¬ 
ative to the Fe sublattice, at low temperatures. Heuce 
a higher temperature is ueeded iu order to reduce the 
Gd momeut to have the same size, albeit with opposite 
directiou, compared to the Fe sublattice. 

It is rewardiug that the agreemeut betweeu theory aud 
experimeut fouud iu Fig. 2 is quite good, aud that the 
three parameters of exchauge iuteractious used iu our 
simulatious explaiu the compeusatiou temperatures of 
the whole rauge of couceutratious of Fig. 2. We also 
uote that fiuer details of the atomic arraugemeut of the 
amorphous structure are very importaut iu achieviug the 
results showu iu Fig. 2, aud this illustrates (as ofteu is the 
case) that atomic arraugemeut (structure) aud maguetic 
properties are iutimately coupled. 


2. All-thermal switching of homogenous samples 

Iu this sectiou, we discuss the dyuamic behavior of 
homogeuous amorphous alloys, usiug the exchauge pa- 
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FIG. 2. (Color online) Curie temperature (Tc) and magnetic 
compensation temperature (Tm) for amorphous Gd-Fe alloys, 
for different concentrations of Gd. The experimental data has 
been taken from Ref. [44]. 



FIG. 3. (Color online) Time evolution of the magnetiza¬ 
tion (M^) for two different concentrations of amorphous Gd- 
Fe alloys under the influence of a thermal heat pulse. The 
solid and dash-dot lines represent the sample concentration 
Gdo. 24 Feo. 76 , while the dash and dot lines are for Gdo.soFeo.ro- 
The magnetization of the two sublattices is plotted separately. 
In the inset we show a typical temperature profile induced 
from the laser fluence, as given by the two temperature model 
of Ref. [31]. 


rameters discussed above, and the magnetic response to 
a femtosecond laser pulse. Only temperature effects from 
the laser pulse were considered, where we adopted a two- 
temperature model, as described in the Appendix C. All 
simulations started with the spin system at room tem¬ 
perature, from which the heat pulse increased the tem¬ 
perature of the spin system in a very short period of time 
(50 fs) to a maximum value, T^ax^ after which the sam¬ 
ple cooled down again. We observed sublattice switching 


for a wide range of concentrations, i.e. 21 at.% to 30 
at.% Gd. As an example we show in Fig. 3 the switching 
behavior of the Gdo. 24 Feo .76 and the Gdo. 30 Feo .70 alloys. 
We find for both concentrations that initially both Gd 
and Fe sublattices demagnetize fast, and that the Fe sub¬ 
lattice reverses its magnetization first, so that for a short 
period of time both Fe and Gd sublattices have parallel 
magnetic moments. Figure 3 shows that after ^3-4 ps 
the all-thermal switching is more or less completed, and 
the reversed sublattice moments relax to their new equi¬ 
librium directions, as the spin system cools down. Alloys 
with different concentrations have slightly different be¬ 
havior, although the main features are independent on 
concentration. All features of Fig. 3 are in agreement 
with the observations of Ref. [1] and also with the simu¬ 
lations and theoretical analysis described in Ref. [9]. 

The initial temperature (300 K) of some of our simu¬ 
lations shown in Fig. 3 is above, or equal to, the com¬ 
pensation temperature (alloys with 21 at.% to 26 at.% 
Gd) while for other simulations the initial temperature 
is below the compensation temperature (alloys with 26 
at.% to 30 at.% Gd). The all-thermal switching behav¬ 
ior was observed for all cases, irrespective whether the 
initial temperature was above or below a compensation 
temperature. The rule-of-thumb proposed for all thermal 
switching in Ref. [45], that it is critical to start with an 
initial temperature below the compensation temperature, 
does not seem to hold in light of the present work. 

In Fig. 4, we show results of the simulated magne¬ 
tization dynamics for various peak temperatures of the 
spin system. Note that if the peak temperature is not 
sufficiently high there is no reversal, the sublattices de¬ 
magnetize to some fraction of their original value, and 
then simply return to the initial configuration, as the spin 
temperature is lowered. This is the case e.g. for Tmax 
= 800 K. If the spin temperature is too high (e.g. the 
simulation with T^ax = 2000 K) both sublattices simply 
demagnetize to a zero moment state, which is stable for 
sufficiently long time in order to make the remagnetiza¬ 
tion completely stochastic in terms of direction of each 
sublattice moment of the final configuration. However, 
for intermediate temperatures (e.g. Tmax = 1000 K) the 
all thermal switching occurs, as is also shown in Fig. 4. 
This shows that it is essential to find the appropriate 
laser fluence with respect to the strength of the exchange 
interactions, for all-thermal switching to occur. 

We end this section with a short note on the effects of 
the damping. We investigated the magnetization dynam¬ 
ics for a wide range of damping parameters, as detailed 
in the Appendix D, and found that the switching behav¬ 
ior reported in Fig. 3 was essentially very dependent on 
the choice of damping parameter. This shows that the 
strength of the intrinsic damping actually determines if 
the all-thermal switching can either occur or not. 
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FIG. 4 . (Color online) Time evolution of the magnetization 
(M^) of amorphous Gd-Fe alloys (Gdo.24Feo.76) for different 
peak temperatures caused by the laser fluence. 


3. All-thermal switching of inhomogeneous samples 

After the initial experimental work of Ref. [1] addi¬ 
tional experimental data was reported, and in particular 
it was argued in Ref. [11] that amorphous Gd-Fe alloys 
may have non-uniform concentration profiles, such that 
some regions are richer in Gd and some are poorer, with 
the opposite trend for the Fe concentration. Interest¬ 
ingly, it was reported that in an all-thermal switching 
experiment of an amorphous Gd-Fe alloy with heteroge¬ 
nous concentration, the Gd magnetic moment reached 
zero before the Fe magnetic moment in the Gd rich re¬ 
gions. After a period of parallel alignment of Fe and Gd 
magnetic moments, the reversal completed with both Fe 
and Gd moments having reversed orientations with re¬ 
spect to their original direction. Hence, Gd rich regions 
exhibited similar behavior as shown here in Fig. 3, albeit 
with the Gd sublattice reaching zero first. 

In order to investigate this experimental result and 
mimic the experimental samples as closely as possible, 
we generated simulation cells with concentration profiles, 
such that, some regions had enhanced (depleted) Gd (Fe) 
concentration with respect to the nominal concentration 
while some other regions had a depletion (enhancement) 
of Gd (Fe) concentration. We considered an amorphous 
alloy with average concentration Gdo. 24 Feo .76 and the Gd 
rich regions had an increase of 6 at.% Gd, while the Gd 
poor regions had a reduction of 6 at.% of the Gd con¬ 
centration. The Fe concentration was modified in the 
same way. A schematic illustration of such a heteroge¬ 
nous sample is shown in Fig. 5, and further details on 
how the inhomogeneous simulation-cells were generated 
can be found in the Appendix E. In these simulations, we 
have in the initial configuration also considered different 


degrees of non-collinearity of the Gd moments. This is 
supported by our first principles calculations, presented 
above, that show an exchange driven non-collinear con¬ 
figuration as the ground state even at T=0 K. From first- 
principles non-collinear theory, we find that the degree of 
non-collinearity is more pronounced for calculations in¬ 
cluding spin-orbit (LS) coupling than without it. More¬ 
over, the effect is more prominent for Gd sublattices (see 
Table III and Appendix A). For example, it is worthy 
to mention here that the current DFT calculations with 
spin-orbit coupling predict a maximal angle deviation of 
Gd atomic magnetic moments with respect to the quan¬ 
tization axis of about 35°. However, the angle deviation 
predicted by DFT without spin-orbit coupling is dras¬ 
tically reduced down to less than the half of the value 
with LS couping. Since the amorphous structure lacks 
inversion symmetry, the Dzyaloshinskii-Moriya (DM) in¬ 
teraction, which favours canted (non-collinear) configu¬ 
rations, is non-zero for both Gd and Fe sublattices. Since 
the spin-orbit interaction is larger for Gd, it is for this 
sublattice we expect a larger effect of the DM interac¬ 
tion. In Table HI we have collected the average and max¬ 
imum angles of the magnetic moments, with respect to 
a common z-axis, for various concentrations. It may be 
seen that exchange effects alone produce a certain degree 
of non-collinearity of the moments, where in particular 
the maximum deviation from collinearity is larger for Gd 
than for Fe. With spin-orbit effects included, the de¬ 
gree of non-collinearity increases, in particular for the 
Gd sublattice. In practice, we have included the de¬ 
gree of non-collinearity, from exchange as well as DM 
interaction, among the Gd atomic moments via a second 
nearest neighbour (SNN) exchange interaction (J^) with 
anti-ferromagnetic character, and as reported below, we 
have followed the simulated magnetization dynamics as 
a function of the degree of non-collinearity of Gd sublat¬ 
tices. 

TABLE III. Degree of non-collinearity predicted by the cur¬ 
rent DFT calculations with and without spin-orbit (LS) cou¬ 
pling. The collected data represent the average (A) and maxi¬ 
mal (M) angle deviation with respect to the z-axis of the mag¬ 
netic moments on Gd and Fe atoms for amorphous Gd^cFei-x 
structures {x = 0.24,0.50,0.76). 


System 


LS coupling 

no LS coupling 


Gd 

Fe 


Gd 

Fe 


A 

M 

A M 

A 

M 

A 

M 

Gdo. 24 Feo .76 

3.2° 

35.1° 

0.6° 1.7° 

0.6 

° 13.6° 

0.4° 

2.2° 

Gdo.5oFeo.5o 

0.6° 

35.6° 

0.6° 4.3° 

0.4 

° 10.3° 

0.3° 

1.4° 

Gdo.76Feo.24 

3.2° 

35.8° 

1.2° 5.7° 

1.2 

° 17.2° 

0.8° 

3.2° 


The results from the simulation cells with inhomoge¬ 
neous concentration are shown in Fig. 6. The figure 
shows results for four different degrees of non-collinearity 
among the Gd atoms in the initial configuration before 
the heat-pulse enters the spin system. The figure also 
shows for each panel the sublattice magnetization of the 
Gd rich-regions and Gd poor-regions. Several things may 
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Fe richGd rich 


Part A Part B 


FIG. 5 . (Color online) Schematic figure showing the inhomo¬ 
geneous samples of Gdo.24Feo.76 plotted in the right side while 
in the the left side is shown a detailed figure of the amorphous 
structure with Fe and Gd atoms in red and blue, respectively. 
Above, Gd rich regions are called Part A and Fe rich regions 
are called Part B. 


be noted from this figure, where the most important fact 
is that the experimental results of Ref. [11] are repro¬ 
duced in these simulations, if a non-collinear configura¬ 
tion of Gd atomic spins are considered within Gd rich- 
regions. Figure 6 shows that for a small degree of non- 
collinear moments of the initial configuration (panels a 
and b in Fig. 6), results in a dynamical behavior that 
is similar to the homogenous results shown in Fig. 3. 
Nonetheless, the data shown in Fig. 3 can be considered 
qualitatively similar to that in Fig. 6 when J 2 is zero, as 
shown by the fact that there is all-thermal switching and 
that the Fe sublattice demagnetizes before the Gd sublat¬ 
tice. Moreover, increasing the degree of non-collinearity 
of the initial configuration causes the Gd moment to de¬ 
magnetize faster, as show in Fig. 6, and for sufficiently 
large values of J 2 it demagnetizes faster than the Fe sub¬ 
lattice. This finding is in agreement with the observation 
reported in Ref. [11]. 

In order to further analyze the results of the heteroge¬ 
nous sample, we show in Fig. 7 the magnetization dy¬ 
namics of the Gd and Fe sublattices, from the different 
regions of the sample, i.e. the Gd poor (Fe rich), Gd 
rich (Fe poor) and regions with an average concentration. 
We can observe that the Gd sublattice switches faster 
(0.42 ps) than the Fe one (1.97 ps) in the Gd-rich regions 
while for the Fe-rich regions the process is reversed. As 
commented above, this finding was already observed ex¬ 
perimentally in Ref. [11] for GdFeCo. The explanation 
proposed by the authors relies on the assumption that 
there are spin currents which transfer torque towards the 
Gd spins in the enriched Gd nanoregions. The theory 
put forth here does not involve explicitly a spin current 
mechanism. Instead, we propose an alternative explana¬ 
tion based on a mechanism driven by a combination of 
the Dzyaloshiskii-Moriya interaction and exchange frus¬ 
tration that produces non-collinearity of the Gd atoms 
belonging to the Gd-rich nanoregion. Using this model, 
the atomistic spin-dynamics simulations show that a non- 
collinear configuration of spins before a heat-pulse enters 



FIG. 6. (Golor online) Magnetization profile (M;^) for Gd 
(blue) and Fe (red) sublattices for inhomogeneous samples. 
The magnetization of Gd rich-regions is shown as solid and 
dash-dot lines (Part A, Gdo. 30 Feo. 70 ) while Gd poor-regions 
are represented by dash and dot lines (Part B, Gdo. 22 Feo. 78 )- 
The strength of the second-nearest neighbour exchange pa¬ 
rameter [J 2 ) for Gd rich-regions is different for any of the 
panels (a, b, c and d) outlined in the figure. For Gd poor- 
regions the J 2 values are considered as 0 mRy. 


the system, explains the faster switching of Gd moments 
in the Gd rich region. We note that the assumption of 
non-collinear moments agrees with the non-collinearity 
predicted by DFT results for the three stoichiometries 
shown in Fig. 1 (see Appendix A). It is important to em¬ 
phasize here the role played by the inhomogeneity of the 












































































FIG. 7. (Color online) Magnetization profile (M^) for Gd and 
Fe sublattices for inhomogeneous concentration profiles in a 
sample with average concentration Gdo. 24 Feo. 76 - The mag¬ 
netization of the Gd rich regions (Part A), Gd poor regions 
(Part B) and sample average is shown in black, green and red 
lines, respectively. The strength of the next nearest neigh¬ 
bour exchange parameter (Jg) for the Gd sublattice was -0.4 
in this simulation. The figure show also similarities with ex¬ 
perimental data reported in Fig. 3c of Ref. [11]. 


sample. Thus, as shown in Fig. 7, only the Gd sublat¬ 
tice in the Gd-rich region switches faster than Fe spins 
while for the sample average, the change of the magne¬ 
tization occurs first for the Fe sublattice. In order to 
measure, detect and use that property for technological 
applications, the experimental techniques are required to 
have at least a nanometer spatial resolution or lower, 
as for example, measuring nanometer-femtosecond spin 
scattering dynamics using X-ray lasers [11]. 

Ultimately, and based on the DFT data collected in 
Table III, we observed that the spin-orbit effects con¬ 
tributes to both sublattices, but are more important in 
Gd sublattices. Moreover, the degree of non-collinearity 
is not evenly distributed over the Gd and Fe atoms as 
shown in Table III by the substantial difference between 
the average and maximal angle deviation of the atomic 
magnetic moment. Consequently, the distribution of the 
DM vectors or SNN exchange interactions is inhomoge¬ 
neous within the Gd-rich and Fe-rich regions. In order to 
mimic and study the effects of the aforementioned inho¬ 
mogeneous distribution, we performed several ASD sim¬ 
ulations with different sets of SNN exchange interactions 
e.g. using values of -0.3 mRy and -0.5 mRy for the Gd 
sublattice distributed randomly in the Gd rich region. 
Interestingly, we observe that an inhomogeneous distri¬ 
bution of the non-collinearity can favour a faster switch¬ 
ing of Gd sublattice if, and only if, the four exchange 
interactions are above a specific threshold (in this set of 


calculations the threshold was -0.3 mRy, as shown for two 
sets of parameters in Fig. 8b). Consequently, an inhomo¬ 
geneous distribution of non-collinear magnetic moments 
can also cause Gd to switch faster than Fe only above a 
minimum value of the degree of non-collinearity. In line 
with the results described in this section, we also stud¬ 
ied the influence of the non-collinearity in the Gd-poor 
region for Gd sublattice and also for Fe sublattice in the 
Gd-rich region. We observe in both cases that the in¬ 
crease of the degree of non-collinearity in this situation 
completely eliminate the switching behavior. We show 
as an example in Fig. 8a the magnetization profile of the 
Gd sublattice in the Gd-poor region for two values of 
the SNN exchange interaction while for the Gd-rich re¬ 
gion the SNN exchange interaction is kept constant as 
-0.4 mRy. The results show that for J^=-0.2 mRy, the 
switching behavior is suppressed for both Fe and Gd sub¬ 
lattices, while if J2 = -0.1 mRy the all thermal switching 
occurs. The conclusion of all these simulations show that 
all-thermal switching is determined by delicate details in 
the concentration profile and the exchange interactions 
of different regions of the sample. 


IV. CONCLUSIONS 

The study of amorphous Gd-Fe alloys represent an out¬ 
standing theoretical challenge because of the lack of crys¬ 
talline periodicity and intrinsic sample inhomogeneities. 
The intricate structural properties clearly determine the 
magnetic ones, as usually is the case, and consequently, 
the magnetization dynamics. We have here been able to 
address the very complicated structural, magnetic and 
dynamical properties of several concentrations of amor¬ 
phous Gd-Fe alloys by using ab initio DFT in conjunc¬ 
tion with an atomistic spin dynamics approach. With the 
aim to assess the validity of this multi-scale approach, we 
compare both structural and magnetic parameters with 
the experimental results and where a comparison can be 
made, we find that they are in a very good agreement 
with observed data. In particular, the Tc and Tm pre¬ 
dicted by our simulations compare quite well with the 
experiment and we are able to explain the increase of the 
compensation temperature with increasing Gd concen¬ 
tration. The explanation mainly resides in the competi¬ 
tion of magnetic sublattices with antiparallel coupling. 

Among the most conspicuous results obtained here, we 
lay emphasis on the crucial role played by the degree of 
homogeneity and non-collinearity of atomic moments in 
the Gd-Fe alloys for the thermally induced magnetization 
switching driven by a femtosecond laser pulse. For ho¬ 
mogeneous samples, the Fe sublattice reverse its magne¬ 
tization before the Gd sublattice for a Gd concentration 
ranging from 21 at.% to 30 at.%. We observe all-thermal 
switching irrespective of whether the initial temperature 
was above or below Tm, which is in clear disagreement 
with previous reported results in literature [45]. In that 
regard, the mechanism proposed for all-thermal switch- 
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Time (ps) 

FIG. 8. (Color online) Magnetization profile (M^) for Gd and 
Fe sublattices for inhomogeneous concentration profiles in a 
sample with average concentration Gdo. 24 Feo. 76 - a) The SNN 
exchange interaction (Jg) in the Gd-poor region (Part B) has 
been chosen to be -0.1 and -0.2 mRy while in the Gd-rich 
region (part A) Jg=-0.4 mRy. b) Two sets of J 2 parameters 
distributed randomly over the Gd-rich region. The values of 
J 2 parameters are listed in the insets of the figure. 


ing put forward in Ref. [45] seems to break down for 
amorphous materials and makes these systems more ver¬ 
satile for spintronic applications since they are less sen¬ 
sitive to the applied initial temperatures. However, for 
inhomogeneous samples, we found the opposite behavior 
with respect to homogeneous case, i.e. the Gd sublat¬ 
tice reaches zero magnetization faster than Fe sublattice, 
at least for the regions with higher Gd concentration. 
Here, we propose a mechanism based on the influence of 
Dzyaloshinskii-Moriya interaction and the exchange frus¬ 
tration that we model by considering a second-neighbour 
exchange interaction between Gd atoms in the Gd-rich 
regions. The microscopic origin of the antisymmetric 
Dzyaloshinskii-Moriya interaction is in general known to 
be coupled to spin-orbit effects and the absence of inver¬ 
sion symmetry, that clearly is present in the amorphous 
Gd-Fe samples. The influence of the damping parameter 
was also considered in this work and we observe that this 
parameter plays a crucial role when dealing with ultra¬ 


fast switching experiments. Thus, the amorphous Gd-Fe 
sample with values of a lower or in the vicinity of 0.02 
undergo a switching process while for higher values of 
the damping, the switching mechanism is totally absent. 
Furthermore, our results point out that a microscopic 
mechanism for all-thermal switching does not need to in¬ 
volve spin current effects. 
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APPENDIX A: ELECTRONIC PROPERTIES 

Total and partial densities of states (DOS) of the Gd- 
Ee system are shown in Eig. 9. The metallic character 
of amorphous Gd-Ee is due to Ee 3(i and Gd bd states, 
which contribute to the DOS at the Eermi level. The 
Gd 4/ states in the LDA+U treatment are localized in 
a narrow energy interval around 8 eV below the Eermi 
level. This is in a rather good agreement with the bind¬ 
ing energy of the occupied 4/ states obtained in XPS 
measurements (9.4 eV) [46]. Although the LDA-bU treat¬ 
ment has been criticized for rare-earths in general [47], 
it is shown from our electronic structure calculations to 
be sufficiently accurate for the purposes of the present 
study. 

Since the electronic structure calculations performed in 
this work have non-collinear spin densities, the spin-up 
or spin-down band picture is no more applicable in our 
treatment. However, the main features of the DOS, i.e. 
the shape of Gd 4/ and Ee 3d states are similar to those 
calculated for ferrimagnetic amorphous Gdo. 33 Eeo .67 al¬ 
loys within DET+LDA theory [48]. Eor both approaches, 
the center of mass of the Ee 3d states is located below the 
Eermi level, while the center of mass of Gd 5d states is 
located above the Eermi level. With the increase of Gd 
concentration, the Ee-Ee bond distance becomes slightly 
shorter, while the Ee-Ee coordination number reduces al¬ 
most 4 times from 7.6 to 2.0. This leads to a narrower 
and less intense Ee 4d states in Gd-rich alloy. We can 
observe, in the middle and lower panels of Eig. 9, that 
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around the Fermi level there is a strong coupling between 
Fe 3d and Gd 5d states. These orbitals are responsible 
for the strong AFM coupling between Fe and Gd atoms. 

The calculated average magnetic moment is 7.40 fiB 
for Gd and 2.38 for Fe atoms. This is in line with 
previous experimental data at T=4.2 K for amorphous 
Gdo. 33 Feo .67 ferrimagnetic alloy [49] and also with theo¬ 
retical values obtained for Gdo. 33 Feo .67 [48]. In Table III, 
we show the average and maximal angles between the z- 
axis and the magnetic moments on Gd and Fe atoms to 
estimate the degree of non-collinearity in the amorphous 
Gda,Fei_a, structures {x = 0.24,0.50,0.76). The results 
exhibit two well-defined features, i.e., Gd magnetic mo¬ 
ments always display angle deviations bigger than Fe 
magnetic moments regardless of the presence of the spin- 
orbit coupling in the calculations, and in addition, spin- 
orbit effects are on the side of an increase of the average 
angle deviation. Some of the Gd magnetic moments are 
much larger than the average value, ranging from 13.6 to 
17.2° without LS coupling and ^ 35° with LS coupling. 
The maximal angle deviation tends to be higher as the 
Gd concentration increases. Thus, the non-collinearity is 
enhanced in Gd-rich structure. Note that our first princi¬ 
ple results contrast the empirical models used to support 
experimental studies on amorphous alloys, which always 
assume the Gd sublattice to be collinear [50]. On the 
other hand, similar as in our results the degree of non- 
collinearity of the Fe sublattice increases with increasing 
the Gd concentration [51]. If the spin-orbit coupling is 
present in amorphous Fe-Gd alloys, as is indicated by 
the current spin-orbit DFT calculations, and due to the 
fact that amorphous structures lack inversion symmetry, 
then these conditions create a suitable environment for 
the DM interaction to be present in these alloys. Even 
though in amorphous materials it is not possible to apply 
straightforwardly the usual symmetry-related arguments 
encompassed by Moriya rules, however it is feasible to 
use the rather general formulas such as the ones derived 
in Refs. [52, 53]. 

Bader analysis shows that in amorphous Gda,Fei_a,, 
there is a charge transfer from Gd to Fe atoms. With the 
increase of Gd concentration from 24 at.% to 76 at.%, 
Fe atoms gain more negative charge, simply due to the 
fact that the probability of finding Gd atoms located in 
nearest neighbour positions around an Fe atom increases. 
At the same time, the average valence electron increases 
by 0.73 and 0.53 for Gd and Fe atoms, respectively. 


APPENDIX B: BINDER CUMULANT 


The fourth order Binder cumulant was introduced in 
Ref. [43] in the context of the finite size scaling the¬ 
ory [54]. For magnetic atoms arranged in a lattice of 
size L, the Binder cumulant is defined by: 


Ul 


<m^ >L 
3 < m? >\ 


(Bl) 



Energy (eV) 


FIG. 9. (Color online) Total and orbital projected densities 
of states calculated for theoretical Gd^Fei-x structure {x = 
0.24,0.50,0.76). The Fermi level is represented by a dashed 
vertical line. 


where m is the order parameter, i.e. the magnetization 
and <> denotes the statistical average taken over sys¬ 
tems at equilibrium and at constant temperature. The 
Binder cumulant allows to locate the critical point and 
the critical exponents in a phase transition. Thus, in the 
thermodynamic limit where the system size of the ferro- 
magnet L ^ 00 and consequently L is bigger than the 
correlation length, the Binder cumulant approximates to 
zero for temperatures higher than Tc while for temper¬ 
atures lower than Tc, Ul —> §• This property of the 
cumulant is very useful for obtaining very good estimates 
of Tc which are not biased by any prerequisites about 
critical exponents. After performing the atomistic spin 
dynamics simulations, we have access to the magnetiza¬ 
tion which is inserted in Eq. (Bl) to obtain as a result the 
Binder cumulant. Then, we plot the cumulant as a func¬ 
tion of the temperature for different sample sizes and Tc 
is estimated from the intersection point of those curves. 
In Eig. 10 we show, as an example, the Binder cumulants 
for two Gd concentrations and the estimated Tc for both 
samples. The Gurie temperatures shown in Eig. 2 have 
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been calculated using the procedure described above. 
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reservoir with eonstant temperature Tiatt- This assump¬ 
tion seems to he quite valid as it was reported in typieal 
pump-probe experiments [55]. Moreover, we also assume 
that the eleetron reservoir is a thermal reservoir mueh 
smaller than the lattiee, hut still larger than the spin sys¬ 
tem. The spin temperature is explieitly passed into the 
stoehastie LLG equation while the eleetron temperature 
ean be expressed in a simple analytieal form as, 

Te = To + (Tp - To) • (1 - (-( 2 ;;^) 

+ (Tp-To)-(l-e(-*/^/)) 

Thus, we redueed the STM model into a simple exponen¬ 
tial funetion [9, 12], that eaptures the essential physies 
of the three temperature model. In Eq. Cl, Tq repre¬ 
sents the initial temperature of the system, Tp is the 
maximum temperature achieved in the simulation and 
Tp is the final temperature of the system. The function 
depends on two time parameters, such as initial time 
and final time r/. The initial time describes the rise-time 
of the temperature to its maximum value and the final 
time represents the relaxation time of the temperature 
from the maximum value to the final temperature of the 
system. Here we used = 50 fs and t/ = 1 ps. In this 
model, the electron temperature is used as spin tempera¬ 
ture. At each time step, the calculated spin temperature 
is passed explicitly into the stochastic field of the LLG 
equation. The values of Tp=800 K, 1000 K and 2000 K 
and To =Tp =300 K are used in the simulations. 


APPENDIX D: DAMPING EFFECTS 


FIG. 10 . (Color online) Calculated Binder cumulants for 
Gdo.24Feo.76 and Gdo.25Feo.75 samples with 1600 atoms per 
unit cell. The size of the samples [x x y x z) was ranged 
from X — y — z— llo A unit cells in steps of 1 . The Curie 
temperature is indicated by the arrows. 


APPENDIX C: DETAILS OF THE 
TWO-TEMPERATURE MODEL 

In order to study the ultrafast demagnetization, the 
three temperature model (STM) was introduced by Beau- 
repaire et al. in 1996 [3]. In the STM model, the electron, 
spin and phonon are thermal reservoirs and coupled to 
each other by coupling constants. It is difficult to define 
an electron temperature for the first femto-seconds of a 
Laser induced pump-probe experiment, but we have here 
for simplicity used the STM model. The analytical ex¬ 
pression of the STM contains three differential equations 
and from that equations, the three temperatures T^ (elec¬ 
tron temperature), Tg (spin temperature) and Tiatt (lat¬ 
tice or phonon temperature) are calculated. In equilib¬ 
rium the temperature of all three thermodynamic reser¬ 
voirs is equilibrated, i.e. Te=Ts=Tiatt- this model 
we assume that the lattiee is an infinitely large thermal 


One of the main parameters in LLG equation is the 
Gilbert damping, a, which is mainly responsible for 
bringing the system into an equilibrium state. It was 
experimentally observed that the damping constant a 
significantly depends on the Gd content and it becomes 
large near to the compensation temperature of the sam¬ 
ples [56]. Though the g-factors of Gd and Fe sublattice 
magnetic moments in our samples are slightly different, 
we used same g-factor for both sublattices, and in our ini¬ 
tial simulations, the static damping parameter was also 
kept equal for both sublattices in the atomistic spin dy¬ 
namics simulations. With these parameters, the evolu¬ 
tion of the magnetization of Gd and Fe moments under 
the influence of an intense femtosecond laser pulse shows 
different precession and reproduce experimental observa¬ 
tions. 

Later on, in a second stage of our simulations, we 
adapted a site-dependent damping parameters in the 
LLG equation and we performed ultra-fast simulations on 
Gdo. 20 Feo. 8 O, Gdo. 24 Feo .76 and Gdo.3oFeo.70 amorphous 
alloys. These results are shown in Figs. 11-13. By fixing 
the damping parameter of Fe species as ape = 0.02, we 
changed the damping parameter of Gd {acd) from 0.02 to 
0.1 in steps of 0.02. If acd = 0.02, the simulations for the 
three sample concentrations shown in Figs. 11-13 predict 
that the Fe sublattice demagnetize faster than the Gd 
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FIG. 11. (Color online) Time dependence of ultrafast mag¬ 
netization (M^) in Gdo. 2 oFeo .80 for different damping param¬ 
eters {a = 0.02, 0.04, 0.08, 0.1) with a fixed electron temper¬ 
ature profile. 


sublattice. Samples with compensation temperatures, 
i.e. with a composition of Gd ranging from 21 at.% to 30 
at.%, shows switching behavior (see Fig. 2). In Fig. 11, 
we found that for Gdo. 20 Feo. 8 O amorphous sample, there 
is no ultrafast switching. This result clearly shows that 
the compensation point is a very important parameter 
in the spin dynamics of Gd-Fe alloy. For higher concen¬ 
trations of Gd, we observed the switching behavior, as 
shown in Figs. 12-13. In our calculations, we found the 
transition metal demagnetize faster than the rare-earth 
element and forms a ferromagnetic-like state for a short 
period of time due to the AFM interaction between Gd 
and Fe atoms. This was already explained in Ref. [1]. 
The idea is that the AFM coupling between Gd and Fe 
atoms favors the spin flipping of Fe atoms when Gd atoms 
are becoming reversed. Thus, the process promotes an 
increase of the net Fe magnetization parallel to the re¬ 
maining Gd magnetization. 

In the case that acd > 0.02, we observed that the Gd 
sublattice moves towards sub-picosecond times but never 
becomes FM to Fe sublattice, as shown in Figs. 11-13. 
The main message of these results is that the damping is 
a very crucial parameter in ultrafast switching process. 


APPENDIX E: GENERATING 
INHOMOGENEOUS SAMPLES 

The inhomogenous sample was constructed from the 
homogeneous unit cell of Gdo. 24 Feo .76 by repeating the 



FIG. 12. (Golor online) Time dependence of ultrafast magne¬ 
tization (M^) in Gdo.24Feo.76 for four different damping pa¬ 
rameters {a = 0.02, 0.04, 0.08, 0.1) with a fixed electron tem¬ 
perature profile. 


unit cell twice in x, y and z direction (12800 atoms), as 
shown in the right panel of Fig. 5. Then, a cube is se¬ 
lected randomly and we reduced the percentage of Fe by 
6 at.%. Thus, the cube turns out to be Gd rich region 
and in the same way the Fe rich regions were created. 
The new supercell consists of an inhomogeneous envi¬ 
ronment and it mimics the original samples of GdFeGo 
experimental results. After that, we study the compen¬ 
sation temperatures for rich and poor areas of Gd. The 
obtained compensation temperatures are similar to ho¬ 
mogenous unit cells. As shown in Fig. 14, we obtained 
a compensation temperature of about 50 K and 350 K 
for Gd poor and Gd rich areas, respectively. The re¬ 
sults shown in Appendix D clearly pinpoint an impossi¬ 
bility of Gd sublattice to switch first than the Fe one. 
Thus, in Sec. IIIG3, we incorporate non-collinearity in 
the sample by introducing an extra AFM exchange value 
to Gd sublattice in the Gd rich areas. Such type of effects 
are observed in the experimental samples. The origin of 
those effects can mainly reside in the concentration of 
Gd-Fe amorphous samples, which modify their structure- 
sensitive properties, such as the magnetic ones, compen¬ 
sation temperatures between the Gd rich and poor re¬ 
gions and also sperrimagnetism found in the Fe sublattice 
of Gd-Fe amorphous alloys for higher Gd-concentration 
while the collinear structure may be expected to exist for 
Gd-poor alloys [51]. 
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FIG. 13. (Color online) Time dependence of ultrafast magne¬ 
tization (M^) in Gdo. 3 oFeo. 7 o for four different damping pa¬ 
rameters {a = 0.02, 0.04, 0.08, 0.1) with a fixed electron tem¬ 
perature profile. 


CQ 

zL 



C 

O 

'-I—* 

CO 

.N 

0 

C 

O) 

CO 



0 100 200 300 400 500 600 700 800 900 


Temperature (K) 


FIG. 14. (Color online) Fe and Gd sublattice magnetizations 
for Fe and Gd rich-regions plotted against the temperature. 
The sublattice magnetization is also shown for regions with 
higher (30 at.%, solid and dot lines) and lower (22 at.%, dash 
and dash-dot lines) Gd concentration. The figure is plotted 
for J2=0 mRy. 
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